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Introduction 


Mathematical modelling of reality usually results in an equation or a set of 
equations. These equations can have a variety of forms, and the nature of the 
solutions to the equations depends on the nature of the equations. 


You should already be familiar with ordinary algebraic equations, whose solutions 
are numbers. For instance, the algebraic equation 

x7-2=0 
has a solution 

x= /2— 1414... 


In Unit 1, you studied equations of a different type called recurrence relations. 
These are equations whose solutions are sequences. For example, we saw that the 
sequence whose nth term u, is defined by 


u, = 2" 
is a solution of the recurrence relation 
Yay =2u, (r= 1,2,3,...) 


The equations we study in this unit are of a different kind again; they are called 
differential equations, Such equations are widely used in mathematical models. You 
will see examples in many units of this course, starting with Units 3 and 4. An 
example of a differential equation is 


dj 

& =x+y () 
You will see that it contains two variables x and y. The problem posed by a 
differential equation is not to find numerical values such as \/2 for the variables in 
it but to find a function relating these variables. One way of describing such a 
function is to give a formula expressing y in terms of x, for example the formula 


y=@-x-1 (2) 


describes one such function. In fact the particular function described by Formula 
(2) happens to be a solution of Equation (1); for if x and y are related by Formula 
(2) then the left-hand side of Equation (1) becomes 

dy _d 


mime Te Hel 


and the right-hand side becomes 

xX+yext(F—x—-l)=e—-1. 
That is, Formula (2) makes the two sides of Equation (1) equal and so the formula 
gives a solution of the differential equation. 


By the end of this unit you will know how this solution was obtained, and you 
will be able to find all the other solutions of the differential equation as well. 


There is an alternative notation for differential equations and their solutions, 
Instead of using two variables, x and y say, this alternative notation uses just one 
variable, x say, together with a symbol for the function the differential equation 
informs us about. For example, if this function were denoted by U, then Equation 
(1) would take the form 


LACM tase) 
dx 


or 
U(x) =x + U(x) 


This is the notation used in 
M101 and MS283, 
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and Formula (2) would be written 
U(x) =e*-—x-1 


or 
Ux —x— 1, 
In this course, however, following the usual practice of applied mathematicians, 


scientists and technologists we shall not normally introduce a symbol for the 
function. 


To translate from one notation to the other, use the following ‘dictionary’ 


this course function notation 
sy U(x) 

dy WU) rg 
dx ae Hie 


Study guide 


There are two main methods for solving differential equations: analytical methods 
which provide the solution function as a formula like Formula (2) above, and 
numerical methods which provide the solution function in the form of a table. In 
addition, there are graphical methods which provide information about the graph 
of the solution function. 


This unit consists of two parts which can be studied independently. One part, 
consisting of Sections 1, 2 and 6, is concerned with numerical and graphical 
methods. It begins with a television programme on some basic ideas about 
differential equations, including a discussion of some graphical techniques; it then 
proceeds, in Section 2, to describe a simple numerical method, Section 6 describes 
the work in your first computer terminal visit, part of which is designed to help 
you learn to use the numerical method described in Section 2. Be sure to plan this 
terminal visit in good time: there may be a deadline after which you cannot use 
the terminal. For up-to-date information about this consult the Course Guide and 
stop press notices, 


The other part of the unit, consisting of Sections 3 and 4, describes some 
analytical methods for solving differential equations. You can do this part of the 
unit either before or after the television programme and computer terminal visit. 
There is a tape commentary which forms part of Section 4. 


Section 5 contains exercises to enable you to revise and consolidate what you have 
learned from the rest of the unit. There is no essentially new material in this 
section. 


The order in which you study the sections of this unit depends on how the 
television broadcast schedule fits in with your study schedule and also on when 
you can arrange your visit to the computer terminal. If you are ready to start 
reading the unit before the television programme is due to be broadcast, begin by 
studying Sections 3 and 4, returning to Sections 1 and 2 after watching the 
programme. Otherwise start by watching the television programme and then study 
the unit in the order it is written. Whichever part of the unit you read first make 
sure you read through the preliminary work at the beginning of Section 1 before 
watching the television programme. Before going on your computer terminal visit 
you must have studied Section 6. This can be done any time after studying 
Sections | and 2. 
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1 Direction fields (Television Section) 


1.1 Preliminary work 


This section describes an approach to differential equations starting from the 
concept of direction field. This concept is explained in the television programme, 
but the following exercise will help you to understand the basic ideas more easily. 


Exercise 1 

On the graph paper provided in Figure 1 draw short line segments 

(i) through (0,0) with slope 1 

(ii) through (3,4) with slope 0 

(iii) through (f, —4) with slope —3 

[Solution on p.43) 

The television programme looks at the practical problem of finding curved mirrors 
which focus a parallel beam of light towards a particular point (see Figure 4). 
Using this problem as an example the programme shows how any direction field 
leads, in a natural way, to a set of curves (called trajectories). The method used to 
find these curves is illustrated by various examples. Exercises 2 and 3 ask you to 
check some of the results used in these examples. 


Exercise 2 
(i) Given that y = 4x? + 5, check that 
dy 
m7* 
(ii) Given that y = 4x? + C, where C is any constant, check that 
dy 
my =x. 


(iii) Given that y = Ce* — x — 1, where C is any constant, check that 


ee 
[Solution on p. 43) 
Exercise 3 
(i) Given that y = (R? — x?) '?, where R is any constant, check that 
d_ x 
ay" 


To keep the calculation under control I suggest the following systematic 
procedure. 


(a) Calculate x from the given formula for y. 


(b) Calculate —~ from the given formula for y. 


(c) Check that the results of (a) and (b) are the same. 


(ii). Given that y= fa —, with C > 0, use a similar systematic procedure to check that 


dy _yt/x?+ 
dx x 
[Solution on p. 43) 


During the programme we shall need to find the orientation of a reflecting surface 
which will deflect a ray of light towards a particular point. To get a formula for 
this orientation we need the law of reflection (see Figure 2). This is an empirical 
law which tells us that the normal (i.e. the perpendicular) to the reflector bisects 
the angle between the incident and reflected rays at the point of reflection. That is 
to say, the angles « and f in Figure 2 are equal. Since « + @ and f + ¢ are also 
equal (both being right angles) it follows that 


0=¢. 
This relation is used as part of a geometrical derivation in the programme. 


reflected ray 


y = (R? — x*)!? is the 
equation of a semi-circle of 
radius R. 


incident ray 


normals, 


point of reflec 


reflector 


Figure 2 
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In order to follow this derivation you also need to know the following two 
geometrical theorems (for how they are used, see Figure 9). 


¥ 
o a eat 8p 


Figure 3(b). If the angles y and 6 are 
equal, then the lengths of sides OQ and 
OP are equal. A triangle such as OOP 
with two equal sides is said to be 
isosceles, 


Figure 3(a). If the lines OQ and PR are 
parallel, then the angles y/ and 0 are 
equal 


Now view the television programme ‘Direction fields and families of curves’. 


° 
1.2 ‘Direction fields and families of curves’ (tv) 3 


In this programme we look at the problem of focusing a parallel beam of light to TVv2 
one point called the focus: we want to formulate a mathematical equation that 
describes a curved reflector which will do this (see Figure 4), 


parallel beam 


{| 


reflector 


Figure 4 Figure 5 


To analyse the problem, we look at one ray of light selected from the beam and 
show that if a small reflecting surface is placed at a point, then there is only one 


orientation of this surface for which the ray is deflected towards the focus (see 
Figure 5). 


By considering other points in the plane, we obtain a whole set of directions. 2 
Every point in the plane (except for the focus) has an associated direction which is You can also think of a 


unique and we can determine it by experiment. Such an association of each point Estevan anion 
with a unique direction is called a direction field (see Figure 6). plane to slopes. 

ee aon hte Gu 

SN SueNe A i tina i 

WS POF ft 

byte Sy OUP ee PP 

Sas, eek tatters eae Se a 

SOE Sine ene) 

Figu 
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At the beginning of the programme, we showed three curves that had the required 
property of focusing a parallel beam of light onto one point, By superimposing 
one of these curves on the direction field (see Figure 7), we immediately notice 
that the slope of the curve at any point coincides with the slope given by the 
direction field. To illustrate how this observation can be used to find equations for 
the curves we first consider a simpler problem. 


We consider the problem of finding equations for the curves that fit the direction 
field given by the formula 


slope at (x,y) = aa 


A diagram of this field is shown in Figure 8. 
4 
y 


4 
i 


es || San 


fon Se ee 


Fo a ewe eee 
7 er I (an i i 
Pape |S ON ON 
(eet a NNN 
f 
i 
ww Sele rf KF 
oy WN NS ee fF 
Se, SN | ae 
Soy SS (a ae te cea d 


Figure 8 


Question 
Sketch on Figure 8 a few curves that fit the direction field. 


This is done by moving your pencil so that it is always travelling in the local 
direction of the field. 


The sketch suggests that the curves fitting the direction field are circles centred at 
the origin, and the computer graphics confirm this; but we would like to be able 
to check this visual judgement by exact calculation. To check that the circles do 
fit, we calculate the slopes of the tangents to the circles. The equation of a circle of 
radius R, centred at the origin, is x? + y? = R? and so expressing y as a function of 
x we obtain (for the part of the circle where y > 0) 


y= (R? — x4)!2, 


Differentiating, we obtain for the slope of the tangent at (x, y) 


SY oh eee ea 
a2 (R? — x2)'2] ~ sy 
So we have checked that the family of circles do fit the direction field whose slope 


at (x,y) is = 
(x, y a 


Now we return to our original problem of focusing a parallel beam of light. The 
example above suggests that it would be helpful to have a formula for the 
direction field. The slope of the direction field at (x, y) is the slope of the mirror 
which reflects the vertical ray of light through the focus O (see Figure 9), By 
considering the triangle PRQ we see that 

PR _PS+SR_y+SR 


OR OR ie 


slope at (x,y) = 


see Exercise 3(i). 


We are using the convention 
that, for example, PR is the 
length of the line from P to R. 
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Since OQRS is a rectangle this becomes 
y+ 09 
= 


slope at (x,y) = (1) 


Now consider the angles. Since OQ and PS are parallel we have (see the theorem 
at the end of Subsection 1.1) 

v=0- 
Also, by the law of reflection we have 

0=¢. 
It follows that 

v=o. 
Thus two angles in triangle OPQ are equal, and so the corresponding sides are 
equal: 

00 =oP 

= Jet ey (by Pythagoras). 


Substituting this into Formula (1) we see that the formula for the direction field is fk | 
| 
i 
i 


[Oa y= / 
+/+) 
slope at -— / 


UJ 
What condition must be satisfied by a curve that fits this direction field? Suppose / 4 
the equation of the curve is @f---------- R 


v=S(x) 


where / is a function we do not know yet. The slope of a tangent to the curve is Figure 9 
dy/dx, the derivative of f(x). If the curve is to fit the direction field this slope must 
match the slope of the direction field, that is, we must have 


dy +/x? +97 


e Beye tes Fs 
An equation like this, connecting the derivative z of some unknown function to 


the variables x and y which are related by it, is called a (first order) differential 
equation. 


There are many different curves that fit the direction field, and each one of them 
corresponds to a different function whose derivative satisfies the differential 
equation. These functions are called solutions of the differential equation. 


Our original reflector problem has now been reduced to the problem of finding 
solutions of Differential Equation (2). At first sight, finding the solutions of a 
differential equation may seem like a formidable task. Fortunately it is often 
possible to find a method which gives all the solutions in terms of a single 
formula. To see how this can be done, we again digress from the reflector problem 
and consider a simpler differential equation 


The solutions of this differential equation can be found by integrating both sides 
with respect to x to give the formula, 


y=h?+C 


where C is an arbitrary constant of integration. Each numerical value for C gives a 
different function relating y to x, and as you checked in Exercise 2 each of these 
functions satisfies (i.¢, is a solution of) the differential equation. Each of these 
solutions corresponds to a different curve fitting the direction field whose slope at 
(x,y) is x (see Figure 10). Figure 10 
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A more complicated example is given by the differential equation 


dy 
ae 


The solutions of this differential equation can be found by the method described in 
Section 4 of this unit. They are described by the formula 

y=Ce*—x—-1. 
Since the process used to obtain this formula is more complicated than a direct 
integration, the arbitrary constant C is this time not just added on to the 
expression for y but is embedded in it in a more complicated way. 


Returning again to the reflector problem, we want to solve the differential 
equation 


dy _ytJ/ety¥ 
dx x , 


This is not one of the standard types considered in this unit, but it happens that 
the equation can be solved by introducing a new variable u, defined by 


u=/x+y 
(The details of the following calculation are not in the television programme and 
are not essential to the course). Squaring both sides of Equation (3) and then 
differentiating with respect to x gives 


du dy 
=2x+ VaR 


2 2 
= 2x + 2y peer) (from the differential equation) 


—(ete + ya/x? +) 
x 


(3) 


2u- 


dx 


2 
= (+4) (by definition of u) 


so that 


= % (from the differential equation) 


It follows by integrating both sides with respect to x that 
u=y+C 
where C is a constant of integration, which must be positive since 
u=/x? + y? > y, Squaring both sides and then using the definition of v gives 
x? + y? =y? + 2Cy + C? 
which simplifies to 


x? = 2Cy 4+ C? 
‘ eo 
ey Sg mas 


This is the formula for the solution of the reflector problem used in the 
programme. 
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Summary of Section 1 
The main points made in the programme are these: 
1. A direction field specifies a slope at each point in the (x,y) plane. 


2. A first order differential equation is an equation connecting the derivative 
dy/dx of some unknown function to the variables x and y which are related by 
that function, 


If y = f(x) is the equation of a curve that smoothly fits the direction field, then 
the function f must satisfy the differential equation obtained by equating dy/dx 
to the expression for the slope of the direction field. Any function (that is, any 
relation giving y in terms of x) which satisfies the differential equation is called 
a solution of the differential equation, There is an infinite family of solutions, 
each one corresponding to a different curve that fits the direction field, 


3. If an exact formula can be found for the solutions of the differential equation 
then this formula will contain an arbitrary constant of integration, called C in 
this programme. Different numerical values for C give different solutions 
and therefore label different curves fitting the direction field. 


4. Whether or not it is possible to find a formula for the solution of a differential 
equation, the direction field is a useful source of information about its 
solution; in particular it can be used to sketch solution curves and so get an 
idea of their qualitative behaviour. 


Exercise 4 

Draw the direction field whose slope is 2 at every point in the plane, Can you deduce the 
family of curves associated with this field? Check your result by differentiation. 

[Solution on p.43] 


Exercise 5 


On p. 8 we checked that the upper half of the circles matched the direction field. Can you 
do the same for the lower half? 
[Solution on p. 43) 


Exercise 6 


What would be the differential equations for the trajectories of the direction fields defined 

by 

(i) Slope at (x,y) = —x/y 

(ii) Exercise 4 

(iii) The direction at any point P (other than the origin) is the same as that of the line OP 
where Q is the origin. 

[Solution on p.43] 


Exercise 7 
For the direction field shown in Figure 11, sketch a few trajectories, including the one that 
passes through the point P. 


Figure 11 2aaetNN Me aN aN 
Ne or as oe, SSE NEN OS EN 
FFE Pe nee ety oe EN NS 
pf ———— 
EL fF 2 oe —-~~>~AN 
[Solution on p.44] f f 4 4 4 a mae 


Exercise 8 
Using the graph paper in Figure 12, draw some slopes belonging to the direction field 
defined by 


slope at (x.y) =1—y 


These curves are often called 
trajectories, 


These functions are sometimes 
called particular solutions, 


This formula is called the 
eneral solution of the 
ifferential equation. 
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including the ones at the points (0,0), (0,4), (0,1), (0, 14), (0,2), (1,0), (1, 1), (1,2). Hence 
sketch enough slopes to give a good idea of how the direction field behaves in the region 
shown, 


Figure 12 


at 


[Solution on p. 44] 


Exercise 9 
Use your direction field from the preceding exercise to sketch some solution curves of the 
differential equation 
dy 
dx 


including the ones passing through the points (0,0), (0,1) and (0,2). 
[Solution on p. 44} 


=l-y 


2 A numerical method 


2.1 A geometrical description of the method 


So far the only systematic method we have for solving differential equations is to 
draw the direction field and sketch its trajectories. In the rest of this unit we 
develop some other methods for solving differential equations. 


For certain types of differential equation there are exact methods of doing this 
which give the solution as a formula such as y = e* — x — 1. Some of these 
methods are considered in Sections 3 and 4; but there are many differential 
equations for which these exact methods do not work and an approximate 
(numerical) method of solution must be used. Most numerical methods yield the 
solution not as a formula but as a table of values, similar to a table of logarithms. 


In this section we shall study the simplest numerical method for differential 
equations. The method is known as Euler’s method, after Leonhard Euler 
(pronounced Oiler), who lived from 1707 to 1783 and was one of the most prolific 
mathematicians of all time (his complete works fill some 60 to 80 volumes). Euler 
devised this method for computing the orbit of the moon; it was used to compute 
lunar tables for the British Admiralty, Euler being voted £300 for his method, 
while the person who did the calculations received £5000. Euler's method is not 
used much now but the more efficient methods that have been developed since are 
based on similar principles. Details of some of these more efficient methods are 
given in the unit on numerical solutions of differential equations. 


Pe Sate C45 ~ 
(a ed a 
= i as 
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A = eg 
~ 
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Le 
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Figure | 
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The idea behind Euler’s method is to approximate the trajectories of the direction 
field associated with the differential equation by segments of straight lines joined 
end to end to form an open polygon as shown in Figure 1. 


Before discussing how to construct this polygonal approximation we must first 
specify which of the infinite family of trajectories we wish to approximate. Let us 
suppose that the trajectory has been specified by saying that it passes through a 
certain point Py with co-ordinates (xo, yo). Such a specification is called an initial 
condition; it provides a starting point from which the polygonal approximation 
can be built up step by step as follows. The first vertex of the approximating 
polygon is taken to be the point Py itself. The first of the line segments forming 
the polygon starts at P, and is taken to have the same slope as the direction field 
at P,; let us call this slope my (see Figure 2). This choice makes the first segment a 
tangent at P, to the trajectory we want, and gives it the right slope at Py; but as 
we proceed along this line segment its slope will deviate more and more from the 
slope of the direction field. If this error is not to become too large, we must soon 
make a change of direction. Let us call the point at which we decide to make this 
change of direction P,: it is the next vertex of our polygon (see Figure 3). At this 
point we begin the next line segment of the polygon, giving it the slope of the 
direction field at P,, which we denote by m,. This new segment is continued to 
some point P, (see Figure 4) at which we once again change direction, the next 
segment having slope m, equal to the slope of direction field at P,. Continuing 
this process we obtain the open polygon shown in Figure 1. 


pe ae ee ion ie on 
is = eer haan ~ 
= ~ ae - ~ ~~ 
al ba ~~ ~ _ a ~ "AS 
aa Pi Ea a end Pp i me 
=a ~ 
os ae a 
~ ~ 
ae Pe ~ a x 
a x x ye 
a = v2 * 
Figure 3 Figure 4 


To follow the trajectory perfectly we would have to adjust the slope of our 
constructed figure continuously to that of the direction field, but in this 
approximation we adjust it only at the points Py, P,,P2,... . 


Exercise 1 


The diagram below shows a direction field whose slope at the point (x,y) is x + y. Using 
the graphical construction just described, draw on the diagram a polygonal approximation 
to the trajectory passing through the origin P,. Choose the points P,,P,,.,.,Ps so that 
their x co-ordinates are 0.2,0.4,..., 1.0, There is no need to make any detailed calculations; 
@ rough sketch is all that is required, 


Po 
a 02 04 6 08 1.0 
wo Se ~-e re 
Ss ~~ e 
Figure 5 


[Solution on p. 44] 
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2.2 A formula for Euler’s method 


If Euler’s method were no more than a graphical construction, it would offer little 
advantage in comparison with sketching solution curves directly from the direction 
field, which is the method we used in Section 1, The point of Euler's method, 
however, is that it can, with advantage, be carried out numerically rather than 
graphically. Instead of a graph showing the polygonal approximation to the 
trajectory, the numerical approach yields a table listing the x and y co-ordinates of 
the vertices Po, P,,..-. By referring to this table we can read off the 
approximation for y corresponding to the value of x we are interested in. 


In order to simplify the calculation of the x co-ordinates we choose the vertices 
Po, P,,... So that their x co-ordinates x9,x,,... are equally spaced, The spacing is 
called the step length and is denoted by h. For example, in Exercise 1 the x co- 
ordinates were 0,0.2,0.4,... and so the step length h was 0.2. In general, suppose 
the trajectory we are trying to approximate has been specified by giving the co- 
ordinates (xo, Vo) of the point Po, If the step length is h then the x co-ordinate of 
P, is given by x, = Xo +h. By adding another step length h we obtain the x co- 
ordinate of P,; x; =X» + 2h, Continuing in this way we obtain the general 
formula, 


X, = Xo + rh 
for the x co-ordinate of the vertex P,. 


The y co-ordinates of the vertices are calculated by setting up a recurrence relation 
from which the y co-ordinate y,,, of the vertex P,,, can be calculated once the y 
co-ordinate y, of the preceding vertex P, is known. To obtain this recurrence 

relation we use the condition formulated in the preceding subsection that the slope 
of the line segment P,P,,, is to equal m,, the slope of the direction field at P, (see 
Figure 6). 


1 
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Figure 6 


¥, 


From the right-angled triangle shown in Figure 6 the condition for the line 
segment P,P,,, to have slope m, is 
length of P,,,Q BS 
length of P,Q i 
Now the length of P,, ,Q is the difference between the y co-ordinates of the 
vertices P,,, and P,, that is y,,, — y,. The length of P,Q is the difference between 


the x co-ordinates of P,,, and P,; it is equal to the step length h. We can 
therefore write 


Ver ~My 


h rare 


Solving for y,,,; we obtain the recurrence relation for the y co-ordinates 
Yer = Y, + hm, 


Since the y co-ordinate of the point P, is known, the y co-ordinates of the 
remaining vertices can be calculated using this recurrence relation. However, 
before we can do this we need a formula for the slope m, of the direction field at 


MST204 2.2 15 


the point (x,,y,). This formula is provided by the differential equation we are 
trying to solve: for example if the differential equation is 


dy 
ant 
then the formula for m, is 
m, =X, + Yer 
For a general differential equation of the type considered in this unit, say 
dy 
ead m(x, y) 


where m(x, y) can stand for any formula involving one or both of the variables x 
and y, the corresponding formula for m, is 


m, = m(x,,Y,). 


To sum up we have: 


Procedure 2.2: Euler’s method 
To apply Euler's method to the differential equation 


dy 
ae = My) 


given the initial condition, y = yy when x = X,, using step length h; 

1. Use the co-ordinates x», yp as a starting point. 

2. Calculate the next y co-ordinate using the recurrence relation 
Jre1 =), + hm, 

where m, = m(x,,),). 

3. Calculate the corresponding x co-ordinate using the formula 
xX, =X trh 

(after increasing the value of r by 1). 

4. If more co-ordinates are required return to Step 2, 


The recurrence relation for Euler’s method can also be derived without reference 
to a graphical construction, By the definition of a derivative, the differential 
equation is equivalent to 
lim 
hoo 
Hence, for sufficiently small h, we have the approximation 


YX + » — yx) _ Preset 


wx +h) = yx) 
h 


This holds for all values of x, in particular for x = x,; so we can replace x by x, in 
the formula and rearrange to obtain 


= m(x,y). 


V(x, + A) = y(x,) + hm(x,, v(x,))- 


With the notation y(x,) = y,, y(x, + h) = y,,, this is the same as the recurrence 
relation in the ‘Euler's method’ box. 


Example 1 

Let us apply Euler’s method numerically to the problem treated graphically in 
Exercise 1. In that exercise we considered the direction field corresponding to the 
differential equation 
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dy 
ie 


and constructed a polygonal approximation to the trajectory passing through the 
origin. The step length was 0.2. 


Since our polygonal approximation started at the origin, we take Py to be the 
origin. So the starting point for our calculation is 


Xp =0,¥9 =0, my =X + Yo =O. 
According to Procedure 2.2, the next y co-ordinate is 
Vi =Yo + 0.2m, = 0 
and the corresponding x co-ordinate is 
xy =X +1x02=02 
which tells us all we need to know about P,: 
x, = 0.2, y, =0, m, =X, +y, =0.2. 
Next we obtain the co-ordinates of P,; 
Ya = +.0.2m, = 004 
and 
X,=Xy+2x02=04 
so we have 
xX, = 0.4, y, = 0.04, m, =X, + y; = 0.44, 
and so on. 


If you have to do more than a very few steps of such a calculation by hand, it is a 
good idea to lay it out as a table, for example: 


Vp m,=%, +), hm, | Yeu = Yt hm, 
0 0 0 0 
0 02 0.04 0.04 

0.04 0.44 


After each value of y,,, has been calculated from the recurrence relation and 
entered in the last column, it is transferred to the y, column in the next row. 


Exercise 2 


Complete the above table and so obtain the coordinates of P,. Check your table by 
comparing it with your graph from Exercise 1. 


[Solution on p. 44) 
Exercise 3 
What changes would be necessary in the column headings of the above table if the equation 


oy 4 yswere replaced by 
dx 


For this new equation, use the same step length, h = 0.2, to calculate an approximation to 
y(1) given that y(0) = 1 (the notation y(1) means ‘the value of y when x = 1’ and (0) means 
‘the value of y when x = 0’). 


[Solution on p.44) 
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2.3 Choosing the step length: accuracy versus cost 


Whenever we use Euler's method we have to choose a value for the step length h; 
in the last subsection this choice was made so as to give a convenient illustration 
of the working of the method. When a method such as Euler's is used for the 
purpose of solving some particular differential equation rather than for illustration, 
however, we have to consider the choice of step length more carefully. Two 
opposing requirements, accuracy and cost, influence this choice, and the decision 
must be a compromise between them. 


The accuracy of Euler’s method usually improves when we take a smaller step, 
size. One way of seeing this is to consider the graphical description of Euler's 
method used in Subsection 2.1. The smaller the step size, the more often we adjust 
the slope of our polygonal approximation to the trajectory, and so the less the 
approximation deviates from the trajectory (see Figure 7). 


=~ as! b= 04 (say) 


- 
re 
~ 
7 
a 
Pe 


Figure 7 exact solution 


The effect can also be seen by calculation. The next table shows approximate 
solutions, using different values of h, for the problem considered in Exercise 3. 


Given that 
dy 


—=y 


dx 
and y(0) = 1, find y(1). 


The table also shows the deviations of these approximations from the exact value 
of y(1), which is e = 2.718282... (you can check that y = e* is a solution of the 
differential equation, that it also satisfies the condition y(0) = 1, and hence that 
y(1) =e! =e), All numbers are given to an accuracy of 6 places of decimals. 


approximation number of 
h to y(1) deviation steps 
+> 
0.1 2.593742 0.124539 10 
0.01 2.704814 0.013468 100 
0.001 2.716924 0.001358 1000 
0.0001 2.718146 0.000136 10000 


You will see that the difference between the approximation to y(1) and the exact 
value decreases roughly in proportion to h. Indeed, it can be shown that by 
making h small enough we can make the error as small as we please, in other 
words, the error approaches the limit zero as h approaches zero. This is true of the 
Euler approximation for any reasonable differential equation, not just the one 
considered in the example. 


A word of caution is necessary here. These remarks about the error are only valid 
if we do the arithmetic to enough decimal places; if we are working with a 
calculator or a computer then the number of decimal places we can use is 
restricted and rounding errors will be introduced. After a certain point any 
increase in accuracy brought about by reducing the size of h will be swamped by 
these rounding errors. 
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Rounding errors are not the only problem. Before concluding that h should always 
be taken very small we must also consider the cost of this additional accuracy. The 
last column of the table illustrates how the number of steps goes up in inverse 
proportion to the step size (the general formula for number of steps is (b — a)/h 
where a and b are the initial and final values of x). Since the error in Euler’s 
method is approximately proportional to step size it follows that for this method a 
10-fold improvement in accuracy is paid for by a 10-fold increase in the number of 
steps required, For the example above, an accuracy of one part in a million could 
require something like a million steps (and each step would have to be calculated 
to a very high accuracy—perhaps one part in 10'*). Even with a high speed 
computer this would not be a sensible way of going about things; instead one 
would use one of the more efficient methods described in the unit on numerical 
solutions of differential equations. Thus Euler’s method is not suitable for high 
accuracy work. Its virtue is rather in its simplicity and its clear illustration of the 
basic principles of the numerical solution of differential equations. 


Exercise 4 
Consider the problem of calculating y(2) (ie. the value of y when x = 2) given that y(1) = 2, 
for the differential equation 
dy 
dx 


using Euler’s method, for various step sizes, (This equation can be solyed exactly to give 
) = }x* + 3, but the purpose of the exercise is to test Euler’s method on a very simple 
example.) 


(i) If the calculation is to contain n steps, what is the step size? 


(ii) Calculate the approximation to y(2) for the step sizes that give (a) I step, (b) 2 steps, 
(c) 4 steps. 


(iii) Calculate y(2) exactly and hence find the errors in cases (a), (b), (c) above. 


(iv) What do the results of (iii) suggest about the dependence of error on h for this 
problem’? 


(v) Assuming your hypothesis in (iv) is correct, how many steps would be required to 
obtain an accuracy of 10~* using Euler's method? 


[Solution on p.44) 


Summary of Section 2 


Given a differential equation 


dy 
i m(x, y) 


and the condition that y(xo) = yp, where x, and yo are specified, then Euler's 
method can be used to find y(b), where b is a given value of x, as follows. 


1, Choose h = (b — Xo)/n where n is the number of steps to be used. 
2. For r=0,1,...,— 1 use the recurrence relation 
Vev1 =I + h(x, y,) 
where x, =X» + rh, Then x, = b and y, ~ y(b). 


The error in the approximation to y(b) is roughly proportional to h, and therefore 
inversely proportional to the cost (ie. number of steps) of the calculation. 


3 Direct integration and separation of 
variables 


3.0 Introduction 


This section and the one that follows can be studied either before or after Sections 
1 and 2. They deal with methods of finding exact solutions, that is solutions 
expressible as a formula, for certain types of differential equation. 
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A feature common to all these exact methods is that at some stage they all involve 
carrying out an integration. We shall see that some types of equation have to be 
rearranged before this integration can be carried out but we begin by considering 
some equations which can be solved by carrying out the integration directly. 


3.1 Direct integration I: general and particular solutions 
Here is one example of a differential equation which can be solved exactly: 


dy 
Frag 


Before proceeding, see if you can find a solution, Write it in the box below, 


y= 


As we have already noted in the introduction to this unit, a solution of a 
differential equation is a function; so we want the box to contain a function of x, 
that is, an expression depending only on x. The condition this expression must 
satisfy is given by the differential equation: it is that differentiation of the 
expression in the box (with respect to x) must give the expression 2x. In other 
words, the box must contain a primitive of the function 2x (a primitive of a 
function f is any function F whose derivative is f), As you probably know one 
primitive of the function 2x is the function x?, and so you could have written x? in 
the box, giving y = x? as a solution of the differential equation. 


This is not the only solution, however; you may recall from an earlier course that 
adding any constant to a primitive gives another primitive. So, for example, the 
function defined by 

y 
is another solution of the differential equation. (You should check this last 
statement.) In fact, any function defined by 


x45 


yex? +C, 
where C is a constant, is a solution of the differential equation, The two solutions 


just considered are the cases C = 0 and C = 5 respectively, but any other 
numerical value for C also gives a solution. 


The method of solution we are using here can be applied to any differential 
equation of the form 

dy 4 
AI”) 
where f is some given function (for example the function 2x in the above example). 
The important thing is that the right-hand side must not involve y. All you have 


to do in such cases is to integrate the function /. The solutions of the differential 
equation are then given by 


y= frosae +¢€ 


where fi (x) dx is one of the primitives of the function f, and C is an arbitrary 
constant. 


The appearance of the constant in this solution is a feature common to all 
differential equations; it expresses the fact that there is an infinite collection of 
functions which satisfy the differential equation, one for each value of the constant 
C. The collection of all the solutions of a differential equation is called the family 
of solutions. 


Just as in the case of recurrence relations, we distinguish between particular and 
general solutions of a differential equation. For the differential equation 

dy/dx = 2x studied above, the function defined by y = x? is a particular solution, 
and the function defined by » = x? + 5 is another. But a formula such as 
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y =x? + C, containing an arbitrary constant, is called the general solution. The 
general solution describes not just one solution but the whole family, from which 
we can pick out any particular solution by giving the arbitrary constant a 
particular value. 


Depending on the problem giving rise to the differential equation, we may require 
either the general or a particular solution. If a particular solution is required, then 
the problem will contain some additional piece of information by means of which 
this solution can be picked out from the family of all solutions. For example, 
suppose our problem is: 


Find a solution of the equation 
dy 

2 

as 

such that (3) = 4. 


The condition y(3) = 4, which means ‘when x = 3, y = 4’, is the additional piece of 
information we need in order to pick out the desired particular solution. The way 
to use this information is first to obtain the general solution of the differential 
equation, which we already know to be 


y=x+e, 


and then use the extra condition y(3) = 4 to give an algebraic equation for C. 
Putting x = 3 in the general solution we obtain 


y=9+C 
and so the condition y(3) = 4 gives 
4=9+C. 


Solving for C gives C = —5 and so the required particular solution is 


xm, 


y 


Exercise 1 
(i) Find the general solution of the differential equation 
toes 
and the particular solution for which y(1) = 5. 
(ii) Find the general solution of the differential equation 
oy 
dx 
and the particular solution for which y(0) = 2. 
(iii) Find the general solution of the differential equation 


mie=2* 


(where a and b are given constants, with b 4 0) and the particular solution for which 
y(0) = 0. 
(Solutions on p. 45) 


3.2 Direct integration II 


All the differential equations we have considered so far in this section have the 
special form 

dy 

oa L(x) 
where f is some given function. Although the derivative of y appears in the 
equation (that is why it is called a differential equation) y itself does not appear, In 
consequence we were able to solve the equation by integrating the function /. We 
can think of this in a slightly different way: in effect, we integrated both sides of 
the equation directly with respect to the same variable x, obtaining 


MST204 2.3 21 


y= fr (x)dx + C N.B. Always integrate both 
sides with respect to the same 
variable, 


where y is a primitive of 2 just as fr (x) dx is a primitive of f, (You might think 


that when integrating both sides of an equation we should also add a constant to 
the left-hand side. However, adding two constants, one to each side of an equation, 
is equivalent to adding their difference, C say, to the right-hand side.) 


Let us now see whether any more general differential equations, containing y as 
well as its derivative, can be solved by integrating both sides of the equation. This equation is equivalent to 
Consider for example, the equation the equation 

dy % = 

2y ih 2x. dx 
; ; discussed in the television 

We know that the right-hand side can be expressed as the derivative of —x? with programme. 
respect to x, The remarkable thing about this equation is that the left-hand side, 
even though it contains y which is related to x in an unknown way, can be 
expressed as the derivative (with respect to x) of a simple expression. Before you 
read on, try to find a simple expression whose derivative (with respect to x) is 
equal to the left-hand side. (Hint: use the chain rule.) 


The required expression is y*; for if we differentiate y*, using the chain rule, we get 
_ 4") dy _ 4 ay 
dx dy “ax ler 


This formula is true whatever functional relation connects y and x, so it does not 
matter that we do not know what this functional relation is. 


The above remarks show that the differential equation can be written in the form 


ay?) _ dx) 
dx dx 


It follows that y* and —x? are both primitives of the same function and must 
therefore be equal or differ by a constant. We can therefore integrate both sides of 
the equation to obtain 


= =-x7+C 


where C is an arbitrary constant. This last formula contains the general solution of 

our differential equation. 

In practice, it is helpful to lay out calculations of this kind as follows. Starting 

from the original differential equation 

dy 

by = 
dx 

we integrate both sides with respect to x, obtaining two primitives which must 


either be equal or differ by a constant. Using indefinite integral notation for these 
primitives, the formula expressing this fact is 


mex 


2 


Remember: there is no need 


See) dx = fi-2max +C, to include a constant on the 
dx left-hand side, one constant is 


Even though we do not know yet how x and y are related, and therefore do not enough. 


know what expressions in x the symbols y and dy/dx stand for, we can use the 
formula for integration by substitution to simplify the left-hand side. In this way 
the equation becomes 


foray = f-2xax +C. 
Writing this formula in the 


Now straightforward integration gives form y? + x? = C shows that 
2 2 the family of solutions can be 
ya —x +, visualized as concentric circles, 
(There is no need to include arbitrary constants when integrating 2y and —2x pee ess aie 


since they have already been taken care of by the constant C.) Since this last programme. 


2 MST204 2.3 


equation shows how the variables x and y are related for each value of C, it 
implicitly describes the family of solutions of our differential equation. However, 
solutions of differential equations are functions and so we must find a formula 
expressing y explicitly as a function of x. This can be done by treating the above 
implicit solution as an algebraic equation and solving it for y to obtain 


y=+/C-x*. 
(The + sign means that for each value of the constant C there are two distinct 
functions satisfying the differential equation: y = ,/C — x? and y = —./C — x*.) 
Notice that in this solution the constant of integration C is no longer added to the 
right-hand side of the equation but has been absorbed under the square root sign. 
Nevertheless, by giving C its various values we still generate the whole family of 
solutions. The above formula is therefore the general solution of the differential 
equation. 


Since the road by which we arrive at such solutions is a fairly long one, the results 
should always be checked by substituting in the original differential equation. In 
the present case, differentiation gives 


Wherever the symbols + or 
appear, the upper sign 
refers to the solution 

so that y = \/C — x? and the lower 
sign to the solution 
y=-/C=¥ 


dy _—— 
we =(+2/C—x) x 


= —2x, 
which agrees with our original differential equation. 


The procedure described above can be summarized as follows. 


Procedure 3.2 
1. This procedure applies to differential equations which can be 
written in the form 


d 
age = Fs) 
mn 

where f and g are given functions. 
2. Integrate both sides with respect to the same variable x 
obtaining 

j (ees frovas $e 

Buy ax ‘i 

Always be sure to include the arbitrary constant C. 


3. The rule for integration by substitution justifies writing this as 


fet vdy= [fends + ¢. 
4. Carry out the two integrations, to obtain a formula relating x 
and y. (This assumes, of course, that the integrations can be done.) 
5. Rearrange the formula in Step 4 to give y in terms of x. 


6. Check that the resulting formula for y satisfies the differential 
equation. If it does then you have found the general solution of 
your differential equation, 


Exercise 2 
In each of the following, ‘solve’ means ‘find the general solution of”. In part (iv) the 


condition (y > 0) ensures that z makes sense; it means that you are only required to find 


solutions for which y is positive for all values of x. 
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(i) Solve (iii) Solve 
dy 1 dy 
2y _ xt eae 
ae i+ydx* 
and check your solution. and check your solution. 
(i) Solve (iv) Solve 
dy 1 ldy 
e—= -—= 0) 
It sian am 


and check your solution. and check your solution. 


[Solution on p. 45] 


3.3 Separation of variables 


This subsection is about a type of differential equation which arises very frequently 
and which can be solved by bringing it to the form considered in the previous 
subsection. One equation of this type is 

d 

9D = x1 +y%), (l) 
dx 
This cannot be integrated directly because the variables x and y appear together 
on the right-hand side, however, we can bring it to the form 


dy, 
av) = LO) 


which can be integrated using Procedure 3.2, by dividing both sides of Equation (1) 
by 1 + y*. This gives 
1 dy 
T+ydx~ 
which is one of the equations you solved in Exercise 2, obtaining the general 
solution 


y = tan (4x? + C), 


This must also be the general solution of Equation (1) for we do not alter the 
information contained in Equation (1) by dividing it on both sides by 1 + y?. 


The procedure we have used here is known as separation of variables, because by 
dividing by 1 + y? we brought the equation to a form where the left-hand side 
depended only on the variable y and its derivative dy/dx, while the right-hand side 
depended only on the variable x: the two variables x and y were separated to the 
two sides of the equation, The procedure can be summarized as follows. 


Procedure 3.3: Separation of variables 


1. This procedure applies to differential equations which can be 
brought to the form 


dy 
ie = FOomy) 
where f and h are known functions. 


2. Divide both sides by A(y), obtaining 


pa igy, 
h(y) dx 
3. Apply Procedure 3.2 to this last equation. 


= f(x). 


A complication which sometimes arises is that h(y) may take the value zero for 
some ‘critical’ value or values of y. Since division by zero is meaningless, a solution 
obtained by the separation of variables method can only be relied on to the extent 
that it avoids such critical values of y. 
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If there are critical values of y (that is, values of y for which h(y) = 0) then there 
will be solutions of the differential equation which cannot be found using the 
separation of variables method. For example, any constant function, with constant 


equal to a critical value of y, will satisfy the equation by making both ws and h(y) 


zero, but to find such a solution using separation of variables would involve 
dividing by h(y) = 0. However, there will also be solutions which avoid the critical 
values of y and they can be found using the separation of variables method, as the 
following example illustrates. 


Example 1 
Consider the equation 
dy 
7 ia (2) 


This has the standard form of Procedure 3.3, with f(x) = 1 and h(y) = y. In this 
case the only ‘critical value’ of y is zero. One solution of the differential equation, 
not given by the separation of variables method, is therefore the constant function 
y =0 (that is, the function which takes the value 0 for all values of x). All the 
other solutions can be found using the separation of variables method, but to 
avoid the critical value 0 we treat solutions for which y > 0 and y <0 separately, 


Dividing both sides of the equation by y (valid for all non-zero values of y) gives 


1 
fiarexee @) 
yi 
where C is an arbitrary constant. Now we want a primitive of 1/y. If y > 0 we can 
use the obvious one, which is log, y, obtaining 
log. y=x+C (y > 0). 
Solving for y then gives 
yoe*e, 
But if y < 0 the primitive log, y will not do since negative numbers do not have 
(real) logarithms. The appropriate primitive of 1/y for negative y is log, (—y). (If 
you have not come across this before you should check that it does have 1/y as its 


derivative, and also note that —y is positive so that it does have a logarithm.) 
Using this primitive in Equation (3) gives 


log.(—y)=x+C  (y <0), 
and solving for y then gives 
y= ete, 


So it seems that the family of solutions of the differential equation requires the 
three formulae 


yet y=-eS y=0 (4) 


for its complete description. However, in Section 4 we shall solve the same 
differential equation by a different method and obtain the general solution as a 
single formula: 


y= Ce (5) 


where C is an arbitrary constant. On the face of it this formula looks quite 
different from Formulae (4) but, as we shall see in Section 4, Formulae (4) and 
Formula (5) describe the same family of solutions. We shall return to this in the 
tape commentary for Section 4 but in the meantime you may like to think about it 
yourself. 


A constant function is one of 
the form y = C where C is a 
number; that is, y takes a 
oneian value C for all values 
of x. 


MST204 23 5 


Exercise 3 


‘Obtain and check the general solutions of the following differential equations. (In some 
cases the allowed values of x and y have been restricted to ensure that you will not be in 
danger of dividing by zero.) 


@ Bee ev 


" dy 
i) xi=y > Oy>0) 


ay Oe 
(iii) ner (y>0) 
[Solution on p. 45} 


Exercise 4 


Obtain and check the general solution of the following differential equations, which are 
used in a later unit. In each case A and B are positive constants. 


. dy 3 d) 

) R=-ay (> 0) Ww) Ae =B-y (B-y<0) 
ti ie py? dy op 

(ii) ae Ay’ (y >0) (v) co ad +y? 


ar dy 
(iii) AT aRay (B-y>0) 


[Solution on p. 46) 


3.4 Partial fractions 

In Unit 3 we shall need the solution of the following differential equation 
dy _ 2 

ax lds 


where a and b are constant. Provided that ay — by? # 0 (that is, y 4 0, y # a/b) we 
can use the separation of variables procedure to bring the equation to the form 
1 dy = 

ay — by? dx 

and integration on both sides then gives 
dy 

oe" 
where C is an arbitrary constant, But to complete the solution we need to evaluate 
the integral on the left-hand side, which is of a type you may not have seen before. 


As always with integrals, we try to transform the integral to one of the types we 
know how to deal with. In this case neither integration by parts nor substitution 
(unless you are extremely clever) is any help; but the following procedure known 
as a partial fraction expansion, solves the problem. 
We notice that the integrand has a denominator which can be factorized: 
1 = il 

ay —by?  (a—by)y" 

The trick is to look for two numbers N, and N, with the property that 


1 Ni NG. 
—— =—1_ 4-2 O and y # a/b). 1 
Gaui rey ies (y# y # a/b), (1) 


Of course it is not obvious a priori that such numbers exist; but if they do then we 
can multiply both sides of the last equation by (a — by)y to obtain the equivalent 
condition 


1=N,y + N3(a — by) (y #0 and y ¥ a/b) 
that is 
1=N,a+(N,—N,b)y (y #0 and y a/b). 


Hint: e**” = e*e” 
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This condition is satisfied for all values of y (except 0 and a/b) if and only if 
1=N,a and 0=N,—N,b. 
This last condition can be solved for N, and N, to give 
N,=1/a and N,=N,b=b/a. 
So provided that a is different from zero our supposition that the numbers N, and 
N, exist is justified and we have, from Equation (1) 
1 bia, Va 
(a—by)y a-by y 
pet ees | 
~ aa — by) © ay 


(y #0 and y # a/b). 


It is a good idea to check such equations by trying a particular value of y; for 
example, substituting y = 1 into the right-hand side gives 
b i ce b+ (a—b) 
a(a—b)" a a(a—b) 
1 

a—b 
which is the same as the expression obtained by substituting y = | into the left- 
hand side. 


Now we can do the integration. To avoid the critical values 0 and a/b it is 
necessary to consider various ranges of y separately. We start by considering the 
case y > 0, a — by > 0: 


dy b 1 
fea ¥ Slams} 
b , 
jSetle 
ajJa—by aly 


= — ioe, (a — by) + Nog, (since y > 0, a — by > 0) 


and so we can finish solving our differential equation. We had already got as far 
as 


dy 
—, = x+C; 
f ay — by? ae 
we can now write 


il 1 
—7 loge (a — by) + logy = x+C. 


To complete the solution we want to solve for y. Multiplying both sides by a and 
using the properties of logarithms gives 

y 
a—by 


log, =a(x + C) 


taking the exponential of both sides and multiplying by a — by gives 
y = (a — by)exp [a(x + C)] 

hence 
y(1 + bexp [a(x + C)]) = aexp [a(x + C)]. 

Finally we have 


__aexp [a(x + C)] 
~ 1+ bexp [a(x + C)] 


a solution which is valid provided y > 0 and a — by > 0, that is, provided 
0 <,y <a/b. Other ranges of y can be treated by the same method. 
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A similar method can be used for any integral where the integrand is the ratio of 
two polynomials. The procedure described below covers only the simplest cases to 
which the method applies. One example of a more complicated case (with a cubic 


in the denominator) is given in Exercise 5(iv). 


Procedure 3.4: Partial fractions 
To express in partial fraction form the function 
ax +p 
ax? + bx +e 


where a, b, ¢, %, are numbers (with a # 0) such that the 
denominator can be factorized into distinct factors. This condition 
will be satisfied if b* — 4ac > 0; for then the quadratic equation 
ax? + bx + ¢ =0 will have two distinct real roots. 


1, Let x,, x2 be the solutions of ax? + bx + ¢ =0. Then the 
function becomes 


eh ies 
a(x — x,)(x — x)" 


2. Suppose that numbers N,, N, exist such that 


ax +B MM 
a(x — x,)(x — x3) 


(x Ax, and x 4 x,). 


1. Sy 
3, Multiply both sides by a(x — x,)(x — x) and collect terms: 
ax + f = aN,(x — x3) + aN,(x —x,) 
=a(N, +N.)x —a(N,x, + Nx) 


4. Since this last equation holds for all x (except x, and x,) we 
must have 


a=aN, +aN, 

B= —ax,N, — ax,N,. 
5. Solve the above pair of simultaneous equations for N, and N, 
6. Check that your result 


(2 N N 
Seat US 
ax? +bx+eo x—x, xX—Xx; 


is correct by substituting one or two particular values of x (not x; 
or xX). 


This procedure does not cover the cases where the quadratic has one or no real 


roots. If it has one real root, say x,, then the substitution x — x, =w will simplify 


the integral. If it has no real roots then a substitution involving arc tan will 


simplify the integral (cf. solution to Exercise 4(v)), But such integrals are beyond 


the scope of this section. 

Exercise 5 

Find partial fraction expansions for 
1 


ar 

ae x+i 

) FaHa 

(ii) Bp : z (where B is a positive constant) 
a 1 

W) Ser 


[Solution on p. 47) 
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Exercise 6 

Using the partial fraction expansions you found in Exercise 5, solve the differential 
equations 

 2ayty — y>o) 

a 1 

@ Qe epeeety 


ii) AM =p y? 
(iii) Az, 7B -y’ (-B<y<B) 
[Solution on p. 47] 


Summary of Section 3 


Differential equations 
The solution of the equation 


#40) 
with fa given function, is 


y= frovas +c 


where C is an arbitrary constant. 


To solve the equation 
a) = S00) 
make y the subject of the formula 
faora = frooax +C. 
To solve the equation 
dy 
dx LAY) th) # 0) 


make y the subject of the formula 


Sion” & frovax +c 


The above are general solutions involving a constant C; to obtain the particular 


solution satisfying a given further condition, use the further condition to evaluate C. 


Integrals 


Lomlansb) ifax+b>0 


f dx _ja 
ax+b 1 F 
qe a= 8) ifax+b<0 


To integrate the function 


ax +p 
ax? + bx +c 


where b? > 4ac, first express it in partial fractions: 


ERA 
X—-X, x—-x, 


where N,, N. are constants and x,, x, are the roots of the quadratic ax? + bx +c. 


You will meet this equation 
again in Unit 4 


Procedure 3.2. 


Procedure 3,3, 


Procedure 3.4. 
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4 Linear equations 


4.0 Introduction 


There is no universal method for obtaining the general solution of a differential 
equation, but the method to be described in this section, the integrating factor 
method, does solve an important class of equations, called linear equations, for 
which the separation of variables method may not work. Like the separation of 
variables method, the integrating factor method depends on rearranging the 
differential equation to a form that can be integrated directly. We begin by 
discussing the type of direct integration appropriate for this. 


4.1 Direct integration using the product rule 


The form of direct integration considered here is based on the formula for 
differentiating a product. The product we are interested in is p(x)y, where p is 
some function to be chosen later, The formula for differentiating p(x)y is 


a x miei ¢ aPlX) 
gy POW] = POT + 


Consequently, if the differential equation we are trying to solve can be written in 
the form 


dy dp(x) _ 
Px) + rey = q(x), (1) 


for some functions p and q, then the differential equation is equivalent to 
4 (p(x)y] = a(x) 

and can therefore be integrated directly, to give 
P(x)y = Java +C. 


Dividing through by p(x) then gives y as a function of x. 


So any differential equation that can be brought to the form of Equation (1) can 
be integrated directly. The next subsection describes how this can be done for 
certain equations. 


Exercise | 


Check that the following differential equations have the form of Equation (1). Determine p 
and q in each case, and hence integrate the equations. 


: dy 
i) nt 


dy F 
(ii) sinx—— + (cos x)y = sinx 
dx 


0 
(ii) eS + Bey =e 


[Solution on p. 48) 


4.2 The integrating factor method 


The method we are going to study here applies to any differential equation of the 
form 
dy 
== 2 
a K(x) + Idy (2) 
where k and | are given functions. An equation having this form is said to be 
linear. 


There is an analogy here with 
linear first order recurrence 
relations, that is, recurrence 
relations of the form 


Yee = ky + ly. 
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To solve Equation (2) we want to find a way of rearranging it to the form 
dy | p(x) _ 
PE +a Y= ae) (1) 


which we saw how to integrate in the previous subsection. (We do not yet know 
what the functions p and q are, but we shall find out during the course of the 
calculation.) As a first step towards bringing Equation (2) to the form of Equation 
(1), we subtract I(x)y from both sides of Equation (2), obtaining 


d 
Se Hoy = Kw. 8) 


To make the first term of this equation match the first term of Equation (1) we 
multiply both sides by p(x), to get 


d 
P(x) — pleylody = plx)ktx). 4) 
The remaining terms of this equation can be made to match Equation (1) by 
choosing the functions p and q so that 
qd 
iets) = =plx)i(x) 6) 


and 
g(x) = p(x)k(x). 


In fact these two conditions are sufficient to determine the two functions p and q. 
We are particularly interested in the function p, because multiplication by p(x) 
takes us from Equation (3), which is not of the directly integrable form, to 
Equation (1), which is. For this reason p(x) is called an integrating factor. 


The condition determining p(x) is Equation (5), a differential equation which can 
be solved by the method of separation of variables. Dividing both sides by p(x) 
and then integrating both sides gives 


fag te = fi-tooas 
that is 
p(x) 
P(x) 
and so one solution is given by 


= fi —I(x)) dx 


log, p(x) = fi —I(x)) dx 


hence 


P(x) = ep| J-onar] (6) 


No arbitrary constant is necessary, nor need we look for solutions with p(x) < 0, 
because we only need one way of bringing our differential equation to the form of 
Equation (1), not all the possible ways. Having determined the integrating factor 
p(x) from Equation (6), we can solve Equation (4) (and hence our original 
differential equation) by direct integration, after which it is easy to find y. The 
following example should clarify the procedure. 


Example 1 
Let us apply the integrating factor method to the equation 


dy _ 5 
moe (7) 


This has the linear form of Equation (2), with k(x) = x and I(x) = 1. The first step 
is to subtract the term containing y, which in this case is y itself, from both sides, 
obtaining 


Equation (7) is discussed in 
the television programme for 
this unit 
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-y=x. (8) 


Now we want the integrating factor p(x). This is best found from Equation (6): 


ple) = exp] f (tea 
=exp| fi-nar] 


=e 


Multiplying both sides of Equation (8) by this integrating factor we obtain 
d 
ere -—e*“y=e 
If the integrating factor has been calculated correctly, the left-hand side of 
Equation (9) should be the derivative of the product P(x)y, that is of e~*y, Using 


the product rule to calculate the derivative of e~*y we find that Equation (9) can 
indeed be written 


(9) 


Lee =e. 
rE [e“*y] =e7*x, 
Direct integration now gives 
e*y= fers dx 
= —xe™* + fersae (integrating by parts) 


== xe" — 67? +6, 


Finally we solve for y, to obtain the general solution of our Differential Equation The solution given in the 
(1): introduction to this unit is a 
particular case of this, with 
y=Ce*—-x—-1. C=1 


The general procedure for solving a first order differential equation by the 
integrating factor method is given as Procedure 4.2 on the next page, 
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Procedure 4.2: The integrating factor method 


1. This procedure applies to differential equations which can be 
written in the form 


dy 
ax 7 KO) + leely 


where k and / are given functions; such equations are said to be 
linear, 


2. Subtract I(x)y from both sides, obtaining 
dy 
ae 7 Oy = Kx), 

3. Define a function p (the integrating factor) by 
pls) = | fir tonya 


(no arbitrary constant is necessary). 
4. Multiply both sides of the equation in Step 2 by p(x), 
obtaining 


dy re 
Plx)o — plx)(x)y = p(x)k(x). 


5. Verify that the left-hand side is the derivative of p(x)y (if not, 
you have made a mistake in calculating the integrating factor) 
and integrate both sides, obtaining 


P(x\y = fi p(x)k(x)dx + Constant. 


Do not forget the arbitrary constant. 
6. Carry out the integration of p(x)k(x). 
7. Solve algebraically for y (that is, divide both sides by p(x)). 


8. Check that the expression you have obtained for y satisfies 
your original differential equation. If it does, you have found the 
general solution of your differential equation. 


4.3 Using the integrating factor method (Tape Subsection) 


The purpose of this subsection is to give you practice in using the integrating \ 
factor method. It consists of three examples accompanied by a tape commentary. 


d 
The first example, the equation 2 = y, has already been solved in Subsection 3.3 


using the separation of variables method; but the method used here is better 
because it covers the cases y > 0, y = 0 and y <0 all together instead of 
separately. 

While working through the examples you will need to refer to the steps in 
Procedure 4.2. 


Start the tape when you are ready, 
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rs] A_worked example 


(2) The integrating factor 


CoS o~- eelren) 


[exe LIC Jax ] 


4 > 
e 


(4) Solve and check 


Gem— , 
d 
Gl & 


sad MST204 2.4 Tape 


P(x) =exp [S- ¢ ddx] 
> = 


¢ 
> 


eee || Sa 
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@ On your own 


Manipulate 
+ identify > 
Steps 1,2 


2a mom = x] 
( Find p(x) 
step 3 


Manipulate. 
steps ae? 


Integrate 
Step 6 


Solution 
Step 7 
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Solution't0 2 my 
dx 


At the end of the tape commentary we asked you to reconcile the general solution 


y= Ae (10) 
obtained by the integrating factor method, with the general solution 
Pers y=-es, y=0 (11) 


obtained using separation of variables. To do this observe that positive values of 
the arbitrary constant A in Solution (10) can equally well be written e“ to give the 
first of Formulae (11), whereas negative values of A can be written as —e to give 
the second of Formulae (11). Setting A = 0 gives the remaining solution y = 0. 


The above example illustrates the fact that it is often possible to find a more 
economical way of writing the general solution of a differential equation simply by 
redefining the arbitrary constant. 


Summary of Secticn 4 

To solve the linear equation 
= K(x) + ey 
dx 


subtract [(x)y from both sides and then multiply both sides by the integrating 


factor 
P(x) = ep| [i-wonac] 


The resulting equation has the form 
dy | dp(x) 
PO +a Y= 9) 
and integrates directly to give the general solution 


PAx)y = favoax +C€. 


5 Revision 
5.1 Choosing the best method 


This unit contains a number of different methods for solving differential equations. 
You will want to know how one decides which to use. 


The first step is to bring the equation to the standard form 


dy 

ros m(x, y). 
Not every differential equation can be brought to this form, for example the 
equation may include second or higher derivatives, but such equations are beyond 
the scope of this unit. 


Next look to see if an exact solution is possible by one of the methods described 
in Sections 3 and 4. In effect, these methods cover two cases only: 


(i) equations soluble by the separation of variables method (Procedure 3.3), 
d 
oe = SOOM), 
ix 
(ii) linear equations using the integrating factor method (Procedure 4.2), 


dy = k(x) + I(x)y, 

dx 
where f and A, or k and J, are known functions, The various cases of direct 
integration can be regarded as special cases of (i) and (ii) (for example, setting 
A(y)'= 1 or I(x) = 0 gives the direct integration considered at the beginning of 
Section 3). 


Equations containing second 
order derivatives are discussed 
in the next unit on differential 
equations. 
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If the exact methods fail, either because m(x, y) has the wrong form or because you 
cannot do the integrations, or even because (as can easily happen) the formula 
given by the exact solution method is too complicated to be useful, you will have 
to fall back on an approximate method. This usually means a numerical method 
and for the time being, therefore, the Euler method. By taking a fairly large step 
size you can carry out the Euler method using a calculator, but if you want 
anything better than a qualitative picture of the solution, you will have to use the 
computer. 


A disadvantage of the Euler method (or any other numerical method) is that it 
only gives particular solutions. If you require information about the general 
solution for a differential equation you cannot solve exactly, you can obtain it by 
plotting the direction field and then sketching enough solution curves to give a 
good picture of their typical behaviour. 


Exercise 1 
For each of the following equations, select the option(s) which give(s) the most appropriate 
method(s) of solution. 


A direct integration 
B separation of variables 
C the integrating factor method 
D_ none of these 
dy 


(i) Fes (vi) 


(ii) a vii) 
Pies ( 


(iii) > = xy (viii) 


(iv) z =x-y (ix) 


(vy) aa 


[Solution on p. 48) 
5.2 Further exercises 


The exercises in this subsection will enable you to revise and consolidate what you 
have learned from the rest of the unit. 


Exercise 2 

For the direction field shown, sketch the trajectory through the point P, 
Wave ON ee AS 
NEVA NEL ARSE Ay AS 
NS Nawks Dy en ek 
NS) ON SOS WOON 
EA ~~N NNN 


[Solution on p. 48} 


Exercise 3 

For the differential equation 
dy 
tee 1—xy 

with (0) = 1, 


(i) write down the recurrence relation for Euler's method with step length 0.1; 
(ii) using step length 0.1 calculate the Euler's method approximation to (0.3). 
{Solution on p. 48) 
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Exercise 4 
If the step size in the above calculation were reduced to 0.05, 


(i) would you expect the error in y(0.3) to increase or decrease, and by what factor 
(roughly)? 


(ii) would the labour of calculation increase or decrease and by about what factor? 
[Solution on p. 48] 


Exercise 5 (Do not attempt this exercise before studying Section 6.) 


If the step size were further reduced to 0.001, write down the options and responses you 
would use to carry out the calculation using the computer program RECREL with print- 
out of the approximations to y(0.1), (0.2) and y(0.3) only. 


[Solution on p. 48] 

Exercise 6 

The general solution of a certain differential equation is 
y=Cx+x?. 


Find the particular solution for which (4) = 8. 
[Solution on p. 48] 


Exercise 7 

Find general solutions of the differential equations 
(i) Fao (x > 0,y > 0) 

(ii) (x > 0) 

(iii) (y>0) 


(iv) 
[Solution on p. 48) 


Exercise 8 
Carry out the integration 
x+4 
[Pears aay 
[Solution on p. 49} 


6 Computer terminal visit for Units 1 and 2 


6.0 Introduction 


See the course guide and possibly also the Stop Press for up-to-date information 
about this visit. Be sure to check whether there is a TMA question requiring use 
of the terminal. 


The purpose of this visit to the computer terminal is twofold; (i) to enable you to 
learn how to use the computer package RECREL to make rapid calculations 
involving recurrence relations and Euler’s method, and (ii) to see some of the 
numerical phenomena described in Units 1 and 2, such as ill-conditioning and the 
balance between accuracy and cost, in an actual computation, 


Before going to the terminal you should have prepared yourself by studying 
Section 6 of Unit | and this section. Section 6 of Unit 1 described in general terms 
how to work the computer packages, and showed you how to turn algebraic 
expressions into valid input expressions for the computer. This section contains a 
description of the package RECREL and some home exercises to help you 
familiarize yourself with it. After completing the home exercises you will be ready 
to start planning your work at the terminal based on the computer exercises in 
Subsection 6.3. It is essential that you prepare these computer exercises carefully at 
home so that you will not waste time at the terminal. 
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6.1 The computer package RECREL 


The computer program RECREL is a package to use when solving first and second 
order recurrence relations. The package can also be used to implement Euler’s 
method to solve first order differential equations. An example of the use of this 
package was given on p. 44 of Unit J. You may need to refer to this to help you 
assimilate the package. 


Command options 

OPTIONS this tells the program to print a list of available options. 

SOLVE this tells the program to run the problem after all the data has 
been input. If some data is missing an error message will be 
printed. 

HELP this tells the program that you are stuck and need advice. 

LIST this tells the program to print the problem and data previously 
input so that you can check which problem you are solving. 

STOP this is the only way to stop the program. 

ROOTS for constant coefficient second order recurrence relations this 


option tells the program to compute and print the two roots of 
the auxiliary equation. 


ANSWER this gives the solution to the computing exercises so that you can 
check that your answer is correct. To the question 


EXERCISE? 
you respond with the number of the exercise you are interested in, 

Problem options 
10 Order of the recurrence relation 
To the question 

ORDER =? 
you respond with either 1 or 2 depending on the order of the recurrence relation 
that you wish to examine. This option is not essential for Euler’s method. 
11 Input recurrence relation 
To the question 

U(R +1) =? 


you respond by typing in the right hand side as a valid input expression (see 
Subsection 6.3 of Unit 1), For example, if you wish to examine the recurrence 
relation 


Ug, = 3u, — 2u,_, + 3r 
you would respond by typing 

3*U(R) — 2*U(R — 1) + 3«R, 
12 Input m(x, y) for Euler’s method 
To the question 

Y'=? 


you respond with a valid input expression for m(x, y). For example if the 
differential equation were 


oy =xy 43 
your reply to Y' = ? would be 
XaY¥ + 3. 


y’ is another way of 


re y 
writing = 
bx 
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Method options 

20 Forward recurrence 

This method uses forward recurrence on either a first or second order recurrence 
relation. The number of terms in the sequence is specified using Option 35. 

21 Backward recurrence 


This method uses backward recurrence on a linear first order recurrence relation 
using Option 32 to specify the starting point. The program will automatically 
adapt the recurrence relation specified by Option 11, For example if you want to 
examine the recurrence relation 


uy, =1—(r + 1)y, 
then you input this recurrence relation using Option 11. The package will 
automatically convert this into the form 
ie Lest 
u,= 
a oo | 


for use with backward recurrence. 


22 Euler’s method 

This uses Euler's method to solve the differential equation 
dy 

Ftd 

where m(x, y) is specified using Option 12. 

Parameter options 

30 Initial condition uy 


This parameter specifies the initial condition u, for first and second order 
recurrence relations. 


31 Second initial condition u, 

This parameter specifies the second initial condition for second order recurrence 
relations. 

32 uy for backward recurrence 

This is used to input the initial condition for backward recurrence with a linear 
first order recurrence relation. N is specified by Option 35. 

33 x» for Euler’s method 

This is used to input the initial value of x for Euler’s method. 

34 y, for Euler’s method 

This is used to input the initial value of y for Euler's method, 

35 The number of terms N in the sequence 

This is used to input the number of terms in the sequence. This must be specified 
for both forward and backward recurrence and for Euler’s method. 

36 Step size h for Euler’s method 

For Euler's method this option specifies the step size. 

Print options 

40 Outline printout 

If you only want selected output you should use this option to specify k such that 
only every kth term in the sequence is printed. 

41 Full printout (default) 

All terms in the sequence will be printed. 


A default option is the one the 
program will use if no other is 
specified 
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42 Print solution only 


Only the last term in the sequence is printed (intermediate terms will not be 
printed), For backward recurrence the value of 1 will be printed, 


6.2 Home exercises 


The following exercises are to enable you to familiarize yourself with RECREL so 
that you will not waste time at the terminal. 


Exercise 1 


(i) It is proposed to put the calculation described in Subsection 2.2 on to the computer; 
that is, to calculate y(1) given that 


dy 
Raxte 


and (0) = 0, with h = 0,2. Write in the table the options you would use and the responses 
to each (where called for) in order to input this problem and print the calculated values of 
(0.2), (0.4),..., v1). 


Information Option Response 


Formula for m(x, y) 
Method 

Intitial value of y 
Initial value of x 
Number of steps 
Step size 

Type of printout 


(ii) Suppose you wanted to increase the accuracy by a factor of 100, so that A is reduced to 
0,002; but that to save printing time you only wanted to print approximations to y(0), 
y(0.1), y(0.2),.... y(1). What changes to the answer to (i) would this entail? 


[Solution on p. 50] 
Exercise 2 
In Sections 3 and 4 of Unit 1 we considered the recurrence relation 
uy) =1—(r+ Du, 
with uo = 1 — e~' = 0.63212056..., which can be used to calculate a sequence of values for 
the integral 
1 
= f eo lytdx, 
° 


Make a table, like the one in Exercise 1, showing the information you must give the 
computer, the options you would use, and the responses to those options, if you wanted to 
calculate the first 25 terms of the sequence, using 


(i) forward recurrence 
(ii) backward recurrence, starting at u,, 
[Solution on p. 50} 


6.3 Computer exercises 


Now here are the exercises for you to try on the computer. Before setting out, 
decide on the order in which you will try them; then decide, for each problem, 
which options you will use and which responses. In this way you will make the 
best use of your very limited time at the terminal. Solutions to these exercises can 
be found from the package using the command option ANSWER 


Note: There may also be a TMA question to be done at the terminal. 


Exercises 3—6 repeat the same numbered exercises from Subsection 6.4 in Unit 1, and 
are included here to refresh your memory. Exercises 7 and 8 you have not met before. 
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Exercise 3 
Use the package to carry out the problem described in Exercise 2 (use both forward and 
backward recurrence), 


Exercise 4 
The recurrence relation 


Wei = (r+ tu, — 1 
with uy =e — 


1,7182818... can be used to evaluate a sequence of values for 
1 
nem f xte!* dx, 
lo 
Use the package to compute the first 25 terms of this sequence. 


Exercise 5 


Check that your mortgage repayments are correct so that you repay the loan in the 
specified time period. (Note that your mortgage may include property and/or life insurance 
which are extras and should be omitted from the calculations.) 


If you do not have a mortgage, work out the monthly payment for a £30,000 mortgage 
borrowed over 25 years if the interest rate is 17°. Use the package to determine the 
amount outstanding at the end of each year, 


Exercise 6 


Use the package to determine the behaviour of the sequences generated by the following 
second order recurrence relations. 


Gi) yay Upaty Wy = 0, 0 1 

(ii) ee = Up — Qty a Wo =0, uy =1 

(iii) 41 =0.9u, — 0.2u,-, + 1500, tip = 3600, wu, = 4000 
(iv) yy = Buy, — 16u,—-, + 4 up = $ u, = 4. 

Exercise 7 


Use the package to carry out the problems described in Exercise 1 (i) and (ii). 


Exercise 8 


This exercise suggests a possible way of treating problems involving differential equations 
which you can not solve exactly. 


(i) Obtain a numerical solution for the reflector equation 


dy _y+J/x?+¥? 
dx x 


which satisfies »(1) = 0, using a convenient step size, say h = 0.2, Take your solution as far 
asx=5, 


(ii) Plot y against x forl <x <5. 


(iii) Try some other step sizes, including h = 0.1, h = 0.05, h = 0.02, and plot a graph to 
show how the computed value of y(5) varies with h. Estimate the exact value of y(5) from 
this graph. 


(iv) (Optional) What would happen if you used a negative h, with the same initial 
condition y(1) = 0? 


The reflector problem was 
treated in the television 
programme, 
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Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


The hypotenuse of the right-angled triangles give the slope of 
the line segments. 


2. In each case we start by differentiating both sides of the 
given formula 


dy d 
® Rabe +5)=JOn+0-% 


fy ence i 
(ii) rier +C)=4(2x)+0=x, 
(ii) We have 
dy d 
(Cr 9 — 1m = 
ak Eee x-1)=C i 
also X+y=x+(Ce*—x-1)=Ce*-1 
dy 
so ae eth 


3. (i) Using the suggested systematic procedure we have 


(a) ae oe — x?)l? = A(R? — x2)" x (—2x) 
-x 
= (R? =x?) 
x —x 
©) Saas 
dy _ x 
(c) a= 5 


(ii) The procedure used in part (i) gives 


(a) +==+0=— 


aN Ee 

(b) yt a +y _2C 2 neg 
~ x 

ie, (oman 

_ 2c 2 4c? 2 4 


From the diagram, the required family of curves appear to be 
straight lines with slope 2, so the family can be described by 
the equation 


y=x+, 
Differentiation gives ¢ = 2+ 0 =2, which checks. 
5. The function giving y on the lower half of the circle is 


y= — faa, 


Hence the slope of the circle at (x,y) is 
== = —4(R? — x2)-"? x (—2x) 

1 x 

=—x (-2x)=—>, 

dy y 
6. In each case the differential equation is found by 
equating ¥ with the slope of the direction field at (x,y) 

d 


Gre 


iD 7e 


x 
y 


Bo Ben 
(ii) oe 


the line OP has slope Z and so the direction field at P also 


has slope 2. Hence the differential equation is 
x 
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Solutions to the exercises in Section 2 
1 


;| 
i a 


The first line segment has slope 0 and ends at the point P, 

= (0.2,0). The next is drawn with a slope equal to that of the 
direction field at P, (the slope can be estimated from the 
diagram of the direction field, there is no need to calculate it) 
and continues until it reaches the point P, at which x = 0.4 
(ie, until it meets the vertical line x = 0.4). The next has 
slope equal to that of the direction field at P, and continues 
until it meets the vertical line x = 0.6, and so on. 


2. The complete table is: 


Yrvn = Yet hn, 


0 

0.04 
0.128 
0.2736 


1 ‘0736 0.21472 | 0.48832 


The co-ordinates of P, are shown in bold type. 


3. The only change to the column headings is that the 
formula for m, is now y, instead of x, + y,. Also, since the 
initial value y(0) is now 1 instead of 0, the first entry in the y, 
column is now 1. The table is therefore: 


0.41472 | 2.48832 


The approximation to (1) is shown in bold type. 


4. (i) The total change in x is 2 — 1 = 1, so with n steps the 
step size is I/n. 


(ii) Euler's formula for the differential equation can be 
written 


VO p41) = YOK) + A, 

so we have 

(a) ifm =1 then h=1 giving 
y2)~ pl) + 1x 1=24+1x1=3 
so for 1 step y(2) ~ 3, 
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(b) ifn =2 then h =} giving 
yd) ~ yl) +4x1=244x1=2 
V2) = y(IQ) +3 x = +x = 3h 
so for 2 steps y(2) ~ 3.25, 


(c) ifn=4 then h =} giving 
y(lt)~2+4% 1= 24 
y(Z) = 2 +4 x = 2s 
(13) = 2y6 +4 x 1b = 28 
y2)= 248 + 4 x 1d = 3h 

so for 4 steps y(2) ~ 3,375. 


(iii) The exact value of y(2) is $x? + 3 with x = 2, that is 
y(2) = 3.5. 


For n= 1 the error is 3.5 — 3 = 0.5. 
For n = 2 the error is 3.5 — 3,25 = 0.25. 
For n = 4 the error is 3,5 — 3.375 = 0.125, 


(iv) It appears that the error is proportional to h, in fact 
the error seems to be equal to 4h, 


(v) For an accuracy of 10~* we would need 4h = 10~*, that 
is h =2* 107*. By (i) the number of steps would therefore 
be 1/h = 5000. 


Solutions to the exercises in Section 3 

1, (i) Integrating the function 6x? gives the general 
solution 

y=wxt+c. 


The required particular solution satisfies y(1) = 5, that is 
2x 1§ + C=5, which implies C = 5 — 2 = 3, The particular 
solution is therefore 


y= 2x9 +3, 
(ii) Integrating the function e~2* gives the general solution 
y= — fe + C. 


The required particular solution satisfies y(0) = 2, that is 
—4e° + C =2, which implies C = 2 +4 = 24. The particular 
solution is therefore 


y= —he-* + 24. 

(iii) Integrating the function asin bx gives the general 
solution 

y= — Fos bx +6. 


The required particular solution satisfies y(0) = 0, that is 


~ Fos (0) + C =0 which implies C = a The particular 


solution is therefore 


a a 
y= = 008 bx + rs 


2. We solve each equation by following the steps listed in 
Procedure 3.2. 


(i) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


fro= fxtax +c 

Step 4: doing the integration gives 

y=be +c 

Step 5: solving for y we obtain the general solution 
y= SEP 43, 


45 


Step 6: check: 


dy 
ax = AGx? + 30)-7? 3x4 = 4 


teas Af at 
: 


(ii) Steps 1 to 3 tell us that we can integrate both sides of the 
equation to obtain 


1 
edy= 
J id ligwere 


Step 4: doing the integrations we obtain 


e” = aretanx + C 
" 


(arctan x means the angle between — 2 


and $ whose tangent 
is x). 
Step 5: the general solution is therefore 
y = log, (arctan x + C), 
Step 6; check: 
dy 


ae xP 3)“ 49, Sx* 


dy 
hence e = 
pal eg 


(iii) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


1 
lnge- xdx +C. 


Step 4: doing the integrations we obtain 
arctan y = 4x? + C, 

Step 5: the general solution is therefore 
y = tan (hx? + C). 

Step 6: check: 


dj 
— = xsec?(4x? +.C) 


dy 
dy xsec? (4x? + C) 
hence —.— = ———_ 7 
1+y?dx 1+ tan? (Ax? + C) 
=x (since sec?) = 1 + tan?()), 


(iv) Steps 1 to 3 tell us that we can integrate both sides of 
the equation to obtain 


1 
~d) = fare 
Se 


Step 4; doing the integrations we obtain 
log,y =x +C (since y > 0). 

Step 5: the general solution is therefore 
y=exp(x+C). 

Step 6: check: 


lay 
ydx exp(x + C) 


x exp(x + C)= 1. 


3. We solve each equation by first applying Procedure 3.3 
and then following the steps listed in Procedure 3.2, 
(i) After dividing by e” the equation becomes 
dy 
yee 
eae 


which has the form described in Step 1 of Procedure 3.2. 


Steps 2 and 3; integrating and using the substitution rule 
gives 


fora = feax +¢. 


Step 4: doing the integrations we obtain 
-ev=e+C. 

Step 5: solving for y gives the general solution 

y= —log,(—e* — C). 

Step 6: check; 

v. 4) x(-e)= =} x (-e)=e", 

(ii) The equation can be brought to the right form for 
Procedure 3.3 by dividing by x giving 

dy 1 


ene 


and we can bring this to the form described in Step | of 
Procedure 3.2 by dividing by y to obtain 
ue 


ydx x 


Steps 2 and 3: integrating and using the substitution rule 
gives 

1 1 

~dy = fea. 

y x 
Step 4: doing the integrations we obtain 
log, y = log. x + C (since x > 0 and y > 0). 
Step 5: solving for y gives the general solution 
y =exp (log, x + C) 


Step 6: check: 
dy. 


(iii) After multiplying by y (or dividing by 1/y) the equation 
becomes 


dy 
I= 


which has the form described in Step J of Procedure 3.2. 


Steps 2 and 3: integrating and using the substitution rule 
gives 


fray = frax +C. 

Step 4: doing the integrations we obtain 

dy? = 4x7 + C. 

Step 5: solving for y gives the general solution 
y= f/x? +2C 


d 
Step 6; check: ve = (x? + 2€)! (x(x? + 2€)"177) 


4. (i) Dividing by y and integrating both sides gives 
fia = fi-avax +C. 
y 


hence log,y = —Ax+C (since y > 0), 


Solving for y gives the general solution 
y = exp(—Ax + C). 


dy 
Check: 7 = —Aexp(—Ax + C) = —Ay. 
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(ii) Dividing by y? and integrating both sides gives 
1 
par fi-Avac+ G 


it 
hence = —-=—Ax+C. 


= 
Solving for y gives the general solution 
1 
=e 
A 

“(ax cF 
(iii) Dividing by B — y and integrating both sides gives 
lr* = fax +c. 


Doing the integrations (using the substitution u = B — y) we 
have 


Check: — Ay? 


- feduex+e 
u 
so —Alog.u=x+C (since u= B—y>0) 
ie. —Alog, (B—y)=x+C. 
The general solution is obtained by solving for y: 
-C€ 
B-y=exp A } 
hence y=B—e>*/e~“/4, 
dj 
Check: A— = e~“/4e-*/4 = By, 
dx 
(iv) ‘The same steps as in (iii) brings us to 
‘A 
= iF du=x+C 
u 
where u = B — y, But since we now have u < 0 the integral 


gives 
—Alog,(—u)=x+C 

ie, —Alog,(—B+y)=x+C. 

The general solution is obtained by solving for y: 


—B+y=exp a4 


A 
fe: y= B+ ere“ Cia 


(er 
Check: AT = ~e"C'4er*'4 = B— y. 


(v) Dividing by B* + y? and integrating both sides gives 


A 
esi facto. 


Doing the integrations (using the substitution y = Btanu if 
necessary) we have 


A y 
Sarctan>=x+C, 
Foon x+ 


Solving for y gives the general solution 


o B(x + C) 
y= Bian ( A } 


dy 
Check: Ap = B? sec? | 


=B ( + tan? Pe+9)) 


B(x + C) 
A 


A 
=B4y? 
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5. We calculate each partial fraction following the steps 
listed in Procedure 3.4, 


(i) Step 1: factorize the denominator to give 
1 1 
+x x(x +1) 
Step 2: assume N, and N, exist such that 
oe ie 
X(x+1) x x+1 
Step 3; multiply both sides by x(x + 1) and collect terms to 
obtain 
L=N,(x +1) + Nx 
=(N,+N)x+N, 
Step 4: it follows that 


Step 5: this is a simultaneous equation with solution N, = 1, 
N, = —1, the partial fraction expansion is therefore 


Al Si 
+x x xt+l 

Step 6: as a check substitute x = 
sides then equal 4. 


(ii) Step 1: the solutions of x? + x 
Hence 


1 into this expansion; both 


Oare 1, —2, 


x+1 x+1 
4x2 DEH) 
Step 2: assume N, and N, exist such that 
x+t x Ny vs Ny ; 
(x—1)(¥+2) x-1 x+2 


Step 3; multiply both sides by (x — 1)(x + 2) and collect 
terms to obtain 


X+1=N,(x + 2)+.N3(x — 1) 
=(N, +N,)x + QN, —N3). 

Step 4: it follows that 

1=N,+N, 

1=2N,— 


Step 5; adding these equations together gives 3N, = 2, that is 
N, =}, substituting this into the first equation gives N, = 4, 
the partial fractions expansion is therefore 


Bee 7 ee 
4x—2 3@—1) | R42) 


Step 6: as check substitute x = 0 into this expansion; both 
sides then equal —4. 


(iii) Step 1: factorize the denominator to give 
1 i 


By” (B-yyB+y) 

Step 2: assume N, and N, exist such that 
1 _ Wi, By 

(B—y)(B+y) Boy" B+y 


Step 3: multiply both sides by (B — y)(B + y) and collect 
terms to obtain 


=N,(B+y)+N,(B-y) 

= (N, —N,)y + (N\B + NB). 
Step 4: it follows that 
O=N,- 
1=N,B+N,B. 
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Step 5: adding B times the first equation to the second gives 
1 =2N,B, that is N, = - substituting this into the first 

1 
equation gives N, = iB the partial fractions expansion is 


therefore 
1 1 1 
By” 2BB—y) * BB+ 
Step 6; as a check substitute y = 0 into this expansion; both 
sides then equal 1/B?. 


(iv) The denominator of this expression is a cubic, but with 
minor modifications Procedure 3.4 can still be used. 
Factorize the denominator to give 


1 1 
x? —1) xe—Nae +1) 
and assume N,, N, and N, exist such that 
SS See 
x(x—1)@+1) x  x- 


Multiply both sides by x(x — 1)(x + 1) and collect terms to 
obtain 


L = Ny (x? — 1) + Na(x? + x) + Ng(x? — x) 
=(N, +N, + N)x? + (Ny —Ny)x — Ny 
Since this must hold for all x (except 0, 1 and —1) we have 
O=N,+N2+Ny 
0=N,-N, 
1=-N,. 
From the third equation, N, = —1, and from the second 
equation N, = N,, These two results, together with the first 


equation, imply that 0 = —1 + 2N,, that is N, = N, = 4. 
The partial fraction expansion is therefore 


Us yer rel aoe 2 
x? 1) x x=)” 2K+1) 


As a check substitute x = 2 into this expansion; both sides 
then equal 4. 


6. (i) After dividing both sides of 


by y? + y, we can integrate and use the substitution rule to 
oul 


vty 
Using the partial fraction expansio >btained in S(i), this can 


be written 

fk-a dy= fate. 

Doing the integrations, we obtain 

log,y —log.(y+1)=x+C (since y > 0). 

By the properties of logarithms, this can be written 

log, al =x+C. 

Taking exponentials of both sides and multiplying by y + 1 
gives 

y=( + Dexp(x + C). 

Solving for y gives the general solution 


__exp(x+C) 
© L=exp(x + Cy) 


exp(x + C) (exp (x + C))? 
l—exp(x+C) (1 —exp(x + C))? 


(ii) Using the partial fraction expansion found in 5(ii) the 
equation can be written 


dy 2 ‘ 1 
dx ~ 3(x —1) * 3(x +2) 
Integrating both sides gives the general solution 


y = flog, (1 — x) + flog, (x + 2)+¢ 
(since x + 2>0 and 1 —x>0). 


dy 
Check: Fi 


Check: 

dy 2 1 
ES aes eet 
dx 3(x-1) 3(%+2) 
Seated 

x +x? 


(iii) After dividing both sides of the equation by B — y?, 
we can integrate both sides to obtain 


A 
erie face, 
Using the partial fraction expansion obtained in Exercise 
S(iv) this can be written 


1 
aj a + 2BB + I 
Carrying out the integrations we obtain 


= fax xe. 


A A 5 
— 55108 (By) + zplog, (B+ y)=x+C 
(since B—y >Oand B+y>0). 

This can be written 


B+y 2B 


ey a +C) 

so (B + y) = (B— y)exp (2B(x + C)/A) 
Bexp (2B(x + C)/A)— B 

hence y= 


“T+ exp (2B(x + C)/A) 


Check: calculation of both ao and B? — y? gives 


4B? exp (2B(x + C)/A) 
(1+ exp (2B(x + ©)/A))™ 


Solution to the exercise in Section 4 


1, (i) This equation can be written in the form 
p(x) 


via +e =H) 


a) 


by choosing p(x) = x and q(x) = 1, for then =LlIt 

follows that the equation can be written 

a(xy) 
dx 

and can therefore be integrated directly to give 


xy=x+C. 


1, 


(ii) Here p(x) = sin x, “a 


= cos x, and q(x) = sin x and 
the given differential equation can be written 
a((sin.x)y) _ 

dx 


It can therefore be integrated directly to give 
(sin x)y = —cosx + C. 


=yty3, 
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(iii) Here p(x) = e* +, = 3e** and g(x) = e* and the 
given differential equation can be written 
de*y) 2, 
ime 
It can therefore be integrated directly to give 


ey = he +0. 


Solutions to the exercises in Section 5 


1. (i) A, (methods B and C could also be used but A is the 
most straightforward), (ii) B or C, (iii) B or C, (iv) C, (v) B, 
(vi) D, (vii) C, (viii) B, (ix) D. 


2 ss 


3. (i) Euler's method gives the recurrence relation 
Veny =¥p +O. x (1 —x,y,) 
where x, =r x 0.1. The method starts from x) = 


Oyo = 1. 
(ii) the calculation is given in the following table 


0.1 
0.089 
0.07622 


The Euler's approximation to y(0.3) is shown in bold-face 
type. 


4. (i) Since the error is roughly proportional to the step 
length, the error would be approximately halved. 


(ii) The number of calculations would be approximately 
doubled. 


5, Option Response 
10 1 
12 1—X*Y 
22 none 
36 0,001 
33 0 
34 1 
35 300 
40 100 


6. Substituting the condition »(4) = 8 into the general 
solution gives 4C + 16 = 8 therefore C = —2 and so the 
required particular solution is 


y= —2x +23, 
7. (i) The standard form of this equation is 
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Using the separation of variables method gives 


1 dy ‘| A 

ae =e (since y > 0) 
1 1 

thus Jee- = frase 

hence =} =-—log.x+C (since x > 0). 
y 


The general solution is therefore 
1 

** log.x—C 

srk 
dx x(log.x—-C x” 

(ii) The standard form of this equation is 
d) 2y 

ae ee 


Check: 


Since this is linear we can use the integrating factor method. 
The integrating factor is 


exp (Res) = exp (2 log, x) = 


Multiplying our original equation by this factor gives 


dy 
eA + 2xy =x? 


which can be written 
ax4y) _ vs 
Sapa 


Direct integration gives xy 
solution is 


= }x* + C and so the general 


(iii) This equation can be solved using separation of 
variables. Since we are only looking for positive solutions we 
can divide by y + y? and integrate to obtain 


fara foree: 


The right-hand side of this equation can be integrated to give 
x + C. The integration on the left-hand side can be carried 
out using partial fractions: let 


1 mths N; 
WO sy 
then 
1=N,(1+y)+ Ny 
=N, +(N, + N2)y. 
It follows that N, = 1 and N, = —1. Hence 


= log, y — log, (1 + y) (since y > 0) 
asi) 
= oes i+y)" 


Equating the two sides of the equation gives 
Delle 
% f + 5] ees 


F y 
Le. ——— ©). 
ie. rarer, exp (x + C), 
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The general solution can therefore be written 
exp (x + C) 
exp(x + C)’ 
2 
Check: 4” = _S*P& + ©) (exp (x + C)) 
dx ~ 1 —exp(x +) ” (1 —exp(x +O) 


(iv) This equation is linear and can be solved using the 
integrating factor method. The integrating factor is 


oo (fi-2a) we 


Writing the equation in the form 
dy 


=yt+y, 


pee 

and multiplying by the integrating factor gives 
dy 

- 2x8 _ o9-25) = _ 9-2 

e & den **y ens, 


which can be written 

d(e~**y) as 
ee 

so enty = 


- fers dx +C. 

Integrating by parts we obtain 

e7 By = xe" — fiertax +c 
= pxe™** + fe* 4. C, 

The general solution is therefore 

y=dr+h+ce*, 


Check: we have 


# =4}+2Ce* 

also 
2y—x=x+442Ce™*—x= 4+ 2Ce* 
and both sides agree, 


8. The integrati 
Since x? + 3x + 


can be carried out using partial fractions, 
= (x + 1)(x + 2) we can write 


x+4 = N, es Ny 
43x42 x+1 x42 
therefore 


X+4=N,(x + 2)+ N3(x +1) 
= (N, + N2)x + (2N, + N2) 


hence 
N,+N,=1 
2N,+N,=4 


The solution of this simultaneous equation is N, 
N2 = —2 and so 


x+4 Pais 2 
xP43x42 x+1 x42 
(Check: if x = —4 then both sides of this equation vanish.) 


We can now carry out the ee 


(ice x+4 «= [Ze 
POR ad +1 se 


= Slog, (—(x + 1)) — 2log, (—(x + 2)) + C, 


where we have used the fact that x < —2 implies (x + 1) and 
(x + 2) are negative, 
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Solutions to home exercises on RECREL in Section 6 


1. (i) The complete table is: 


Information Option Response 
Formula for m(x, y) X+Y 
Method (Euler) none required 
Tnitial value of y O 
Initial value of x ) 
Number of steps 5 
Step size 02 


Type of print-out none required 


(ii) The response to Option 36 becomes 0.002; Option 4 | is 
replaced by Option 40, and the response to it is 50 because 
we only want every 50th calculated value (since 50 x 0.002 

= 0.1, the interval between successive printed values of x). 


2. (i) A suitable table is 


Information Option Response 
Order of recurrence relation 10 1 
Recurrence relation is 1 1 —(R + 1)*U(R) 
Method (forward recurrence) 20 none required 
Initial value uy 30 0.63212056 
Number of terms 35 25 
Type of print-out 41 none required 


(ii) Part (i) can be left in the machine but the following new 
options must be inserted. 


Information Option Response 
Method (backward recurrence) 21 none required 
Value of uy 32 any number 
between 0 
and 1 say 
Number of terms 35 35 


Note: Inserting Option 21 automatically cancels Option 20 
given in part (i), also, the response 35 to Option 35 
automatically cancels the previous response of 25 to Option 
36 given in part (i). 


